LISA observations of supermassive black holes: parameter estimation using full 

post-Newtonian inspiral waveforms 
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We study parameter estimation of supermassive black hole binary systems in the final stage of 
inspiral using the full post-Newtonian gravitational waveforms. We restrict our analysis to systems 
in circular orbit with negligible spins, in the mass range 10* Mq — 10^ Mq, and compare the results 
with those arising from the commonly used restricted post-Newtonian approximation. The conclu- 
sions of this work are particularly important with regard to the astrophysical reach of future LISA 
measurements. Our analysis clearly shows that modeling the inspiral with the full post-Newtonian 
waveform, not only extends the reach to higher mass systems, but also improves in general the 
parameter estimation. In particular, there are remarkable improvements in angular resolution and 
distance measurement for systems with a total mass higher than 5 x 10^ Mq, as well as a large 
improvement in the mass determination. 

PACS numbers: 04.25.Nx, 04.80.Nn, 95.55.Ym, 97.60.Lf 



I. INTRODUCTION 

Supermassive black hole binary systems, in the mass 
range lO^M© — lO^M©, will be detectable by the Laser 
Interferometer Space Antenna (LISA) [l], [3, [1] through- 
out the entire Universe. Observations of gravitational 
waves from this class of sources are among its high- 
est priority targets. By measuring these gravitational 
waves we will have detailed information regarding gen- 
eral relativity itself and the behavior of space-time 
[1, H, H, 0, [1, 3 ; precision measurements of the Universe 
as a whole Ell, [H, [3; the formation and growth 
of massive black holes in galaxy evolution [ll|, [3, [l^l ; 
and black hole astrophysics [l^, [l3| . Detection of grav- 
itational waves will provide complementary information 
to conventional astronomy. 

Supermassive black hole binaries are long lived sources 
in the LISA band. The whole coalescence of a compact 
binary system is commonly divided into three phases: 
the adiabatic inspiral, the merger and the ringdown. 
Most of the signal-to-noise ratio (SNR) accumulates dur- 
ing the last days prior to coalescence and during the 
merger phase, but one critically relies on long integration 
times to disentangle the source parameters, in particular 
to resolve the source position in the sky and measure its 
luminosity distance [18| . This is due to the motion of 
LISA around the Sun that breaks the degeneracy in the 
parameters. Thus, in this paper we restrict our attention 
to the inspiral phase since this is the most interesting for 
parameter estimation. 

Because of the complexity of the problem, most anal- 
ysis carried out so far to address how accurately LISA 
can measure the source parameters and the implica- 
tions for astronomy and cosmology, have considered 
only the inspiral restricted post-Newtonian waveforms 
0,llll,[ll,[l3, 18, 19, 20, 21, 22, 23], where all amplitude 
corrections are discarded and only post-Newtonian (PN) 



contributions to the phase are taken into account. Thus, 
the restricted-PN waveform consists of just the domi- 
nant harmonic at twice the orbital phase. Only in a few 
cases, it has been reported the importance of including 
higher order PN terms to the amplitude and the phase 
0, m, m, [13 • In the context of ground based detec- 
tors, it was found that the consequences of arnplitude 
correction in the templates are considerable [H, [2^, [s^l . 

The main purpose of the present analysis is to inves- 
tigate in detail the impact of full versus restricted PN 
waveforms for parameter estimation, by exploring a vast 
parameter space, specifically in the context of LISA, ex- 
tending previous results. The waveform we use is de- 
scribed by 11 parameters. Therefore extensive large- 
scale, CPU-intensive Monte Carlo simulations have been 
required for such an exhaustive study. 

The conclusions of this work are particularly impor- 
tant with regard to the astrophysical reach of future 
LISA measurements. Our analysis clearly shows that 
modeling the inspiral with the full post-Newtonian wave- 
forms, as compared to the restricted-PN ones, not only 
extends the reach to higher mass systems, up to 10^ Mq, 
as discussed by Arun et al. [2^, but also improves the 
parameter estimation. Improvements in angular resolu- 
tion and distance measurement are remarkable for sys- 
tems with a total mass higher than 5 x 10^ Mq. These 
results are in agreement with those recently found in . 

The organization of this paper is as follows. In Sec- 
tionini we discuss the LISA detector output, and we pro- 
vide the total LISA noise curve employed in our analy- 
sis. In Sec.|lll]we describe the gravitational wave signals 
from binary systems. Sec. II VI reviews the basic concepts 
of signal parameter estimation in matched filtering. In 
Sec.|V]we spell out the assumptions about the waveform 
and the observations on which our analysis is based, we 
present a detailed description of the waveform model and 
all key steps to compute the SNR and the estimation of 
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the parameter errors. Sec. IVII presents the results of 
our investigations, where we compare the impact on pa- 
rameter estimation of supermassive black hole binaries 
using full post-Newtonian inspiral waveforms versus the 
restricted ones. The results are presented for a given 
source location, exploring a vast parameter space, we 
also study how source location errors in terms of ad- 
vanced warning times, and the errors dependency with 
redshift. Finally, Sec. IVIII concludes with a summary of 
the main results of this paper and present pointers to 
future work. 



II. LISA RESPONSE TO GRAVITATIONAL 
WAVES AND DETECTOR NOISE 

LISA consists of three spacecrafts arranged in an equi- 
lateral triangle orbiting the Sun. The arms of the trian- 
gle are approximately L = 5 x 10^ km in length, and the 
triangle is inclined at an angle of 60° to the ecliptic. The 
entire triangular configuration spins as the antenna or- 
bits the Sun, rotating once during a single orbit. A grav- 
itational wave interacting with the configuration causes 
the length of the three arms to oscillate. 

For the supermassive black hole binary inspirals con- 
sidered in this paper, most of the SNR accumulates at 
frequencies / < 10 mHz, so it is adequate to use the low- 
frequency approximation to the LISA response function 
derived by Cutler p^ . In this approximation, LISA can 
be regarded as two independent gravitational wave de- 
tectors with 90° arms and rotated 45° with respect to 
one another. 

The strain h{t) produced at the output of LISA 
Michelson interferometer by a gravitational wave (GW) 
signal characterized by two polarization states (t) and 
hx{t) is 



(2.1) 



where and i^^*' are the time-dcpcndcnt antenna pat- 
tern functions, the factor ^/3/2 comes from the 60° open- 
ing angle of the LISA arms, and the i =1,11 labels the 
two independent Michelson outputs. The response func- 
tions F^\t) and F^^^t) depend on the direction and 
orientation of the source in the sky and they vary with 
time because during the observation, the interferometer 
changes orientation with respect to the source. Wc refer 
the reader to [ll] for further discussions and details. 

The total noise that affects the observation of radia- 
tion emitted by binary systems is given by the superpo- 
sition of instrumental sources, S'^''''(/), and astrophysi- 
cal foregrounds of unresolved radiation, S'™"'(/), the so- 
called confusion noise. The total noise spectral density 
Snif) is therefore the sum of these two components 



Snif) = s'rif) + sr'if)- 



(2.2) 



A good fit for the instrumental noise, for / < 5 mHz, 
is given by [3l|, [111 

ST'if) = 6.12 X 10-^1 f-^ + 1.06 X lO-''^ 

+6.12 X 10"^^ f Hz"^ (2.3) 

where / is in Hz. This is derived from the online sensi- 
tivity curve generator (ssj . which is based on the noise 
budgets specified in [l|. 

For the confusion noise we adopt the same analyti- 
cal approximations given in |3l| . considering only noise 
from short-period galactic 5','^™° and extragalactic bina- 
ries 5*^^° (due to white dwarfs binaries) , assuming they 
are all unresolvable, i.e.: the worst case, and we ignore 
the effects of captures of compact objects. That is, we 
estimate the total effective noise density as 

5,f(/) = min{[5;r(/) + Sr°(/)] exp(AcT-idA^/d/), 

^r(/) + 5r°(/) + sr°(/)} (2.4) 



where we take kT ^ = 1.5 yr ^, 

^ = 2 X 10-3 /-"/3 Hz-i 
d/ 



Sn'^^if) = 1.4 X 10-** / 



-7/3 



Hz" 



and 



III. 



S'^™(/) ^ 2.8 X 10-*^ Hz" 



(2.5) 



(2.6) 



(2.7) 



BINARY BLACK HOLE COALESCENCE 
WAVEFORMS 



The noise contributions in each data stream I and II are 
by definition the same. 



The coalescence of binary black holes is commonly di- 
vided into three successive epochs in the time domain: 
inspiral, merger and ringdown. During the inspiral the 
distance between the black holes diminishes and the or- 
bital frequency sweeps up. The waveforms are well mod- 
eled using the post-Newtonian approximation to general 
relativity. Eventually the post-Newtonian description of 
the orbit breaks down, and the black holes cannot be 
treated as point particles any more. What is more, it is 
expected that they will reach the innermost stable circu- 
lar orbit (ISCO), at which the gradual inspiral ends and 
the black holes begin to plunge together to form a single 
black hole. This is referred as the merger phase. At the 
end, the final black hole will gradually settle down into 
a Kerr black hole. 

In recent years, a series of breakthroughs has occurred 
in numerical simulations of binary black hole systems 
[3^ . Issl . |36| . Long-term evolutions of inspiralling black 
holes that last for several orbits have been obtained with 
several independent codes [13, [H, M, \M, \M, [H, [3, 
[45I [46t , and accurate GW signals have been computed, 
including the merger and ringdown phases. Still, the 
post-Newtonian approximation to general relativity is 
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the best available method for calculating the vast ma- 
iority of the GW signal cycles observed by LISA. See 
[47[ for a review and extensive references. 

For non-spinning black holes, the best PN waveforms 
currently available have been calculated at 2.5PN order 
in amplitude and 3.5PN order in phase [H, |4^, \E^, [sil . 
[s^ . [531 ■ These waveforms in the two polarizations h+ 
and hx take the general form 

+vy'Hi'/^^+v'H^^^,+v'/'H^l/^^} , (3.1) 

where we have set G = c = 1, as we will do throughout 
this paper, v = (Muj)'^^^, uj is the orbital frequency, Dl 
is the luminosity distance to the source, and M and 77 are 
the observed total mass and the symmetric mass ratio 
respectively, defined in Sec. [V] The explicit expressions 
for H^l^^ , m 0, . . . , 5 can be found in 0, |^ ■ They 
include contributions from several harmonics of the bi- 
nary's orbital motion. 

For black holes with significant spins, the state of 
the art is somehow less advanced and the correspond- 
ing waveforms have been calculated through 2.5PN order 

Equation (|3.ip corresponds to the so-called full wave- 
form (FWF). Given its complexity, together with the fact 
that the second harmonic contributes most strongly to 
the waveform over most of the inspiral phase, it is com- 
mon to make some simplifications and work only with 
the restricted waveform (RWF), in which one neglects all 
amplitude terms except the Newtonian quadrupole one, 
but keeping the phase to some specific PN order, i.e., 



keeping only i?,^" ^ and throwing out the rest H 
m > 0. 

It is the goal of this paper to revisit the problem of pa- 
rameter estimation for supermassive black hole binaries 
studying the improvement in error estimation by using 
FWF and compare with the previous results obtained 
with the RWF. 



Ani/2) 



for 



IV. REVIEW OF SIGNAL ANALYSIS AND 
PARAMETER ESTIMATION 

In this section we briefly review the basic concepts and 
formulas of signal parameter estimation relevant to the 
goal of this paper; we refer the reader to [13] for a more 
detailed analysis. 

The signal s'^*'(i) as measured by the detector i is 
a superposition of noise n^''\t) and gravitational waves 
/i«(t;A) 



gW(f) = /iW(i;A)-Hn(')(i) 



(4.1) 



where A represents a vector of the unknown parameters 
(location, masses, spins, etc) that characterize the actual 



waveform and that one wishes to estimate from the data 
stream. 

For sake of simplicity we shall made the standard as- 
sumptions that the noise n'^^\t) has zero mean and it 
is stationary and Gaussian. Within this approximation, 
the Fourier components of the noise are statistically de- 
scribed by 
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E[n{f)n*{n]^M-nSn{f) 



(4.2) 



where E\\ denotes the expectation value with respect to 
an ensemble of noise realization, the * superscript de- 
notes complex conjugate, 5'„(/) is the one sided noise 
power spectral density, and tildes denote Fourier trans- 
forms according to the convention 



i(/) = / e'^'^f'xit) dt . (4.3) 



With a given noise spectral density for the detector, 
one defines the "inner product" between any two signals 
g{t) and h{t) by 



{g\h)^2 



~9*{f)Hf) + ~9{f)h*{.f) 
Snif) 



d/. 



(4.4) 



With this definition, the probability of the noise to have 
a realization hq is just 



p{n = no) oc e-("''l"°)/2 , 



(4.5) 



The optimal signal-to-noise ratio (SNR) p, achievable 
with linear methods {e.g., matched filtering the data) is 
given by the standard expression 



(4.6) 



In the limit of large SNR, which is typically the case 
for LISA observations of supermassive black hole binary 
systems, the probability that the gravitational wave sig- 
nal h{t; A) is characterized by a given set of values of 
the source parameters A = {A'"'} is given by a Gaussian 
probability of the form [sst 



p(A|/i) =p("'(A)cxp 



1 



(4.7) 



where AA*' is the difference between the true value of the 
parameter and the best-fit parameter in the presence of 
some realization of the noise, p^") (A) represents the dis- 
tribution of prior information (a normalization constant) 
and r^fe is the so-called Fisher information matrix de- 
fined by 



(4.8) 



dih*{f)dfh{f)+d,h{mh*{,f) 

Snif) 



df: 
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where di = 

The inverse of the Fisher matrix, known as the 
variance-covariance matrix, gives us the accuracy with 
which we expect to measure the parameters A*^ 

j.jk ^ (r-i)jfc ^ (AA^AA*^) . (4.9) 

Here the angle brackets denote an average over the prob- 
abiUty distribution function in Eq. (|4.7p . The root- 
mean-square error ak in the estimation of the parameters 
A*"' can then be calculated, in the limit of large SNR, by 
taking the square root of the diagonal elements of the 
variance-covariance matrix. 



a,. = ((AA'=f)V2 = 




(4.10) 



and the correlation coefficients c^'^ between two param- 
eters A^ and A*' are given by 

,,, (AA^AA*^-) EJ'^ 
^3L^\ i = (4.11) 

Returning again to the two detector case, we have that 
the largest value of the SNR is 

r = ^{p'Y + {p'')\ (4.12) 

{p'^"f = {h'^"\h''''), (4.13) 

and we can write a total Fisher matrix as the sum of the 
individual Fisher matrices for each detector 

r-'-r(,+r(/. (4.i4) 

In this paper wc will use equation (j4.9p . along with 
the FWF and RWF models for the binary black hole 
coalescence, and the LISA noise spectrum discussed in 
Sec. ini in order to estimate how well LISA will be able 
to measure the source parameters. We refer the reader 
to [5^, [B^I for a detailed discussion about the conditions 
required for the Fisher-matrix formalism to be applica- 
ble. 



V. OBSERVATION OF SUPERMASSIVE 
BLACK HOLE BINARY SYSTEMS 

We consider observations of supermassive black hole 
binary systems of masses mi and TO2 at luminosity dis- 
tance Dl. For later convenience we define the following 
mass parameters: total mass M = mi -\- m2, reduced 
mass fi = mim2/M, the symmetric mass ratio rj = fi/M, 
and M = ^^/^Af^/^ = Mrf/^ the chirp mass. We focus 
on binary systems in the mass range IQ^Mq — IQ^Mq 
and wc do not consider here the case of binaries with an 
extreme mass ratio, e.g, a black hole of lOM© orbiting a 
supermassive one. The reason is that some assumptions 
about the waveform that we will be considering would 
be rather unrealistic for such astrophysical system. 



Before presenting the results, in this section we spell 
out the assumptions about the waveform and the obser- 
vations on which our analysis is based, together with a 
detailed description of the post-Newtonian inspiral wave- 
form we use. 



A. Assumptions 

The waveform model we consider is based on the fol- 
lowing assumptions: 

• As signal we consider only the inspiral phase of 
the whole coalescence, neglecting all information 
coming from the merger and the ringdown phases. 
We terminate the inspiral when the binary's mem- 
bers are separated by a distance 6M; this very 
roughly corresponds to the point at which the post- 
Newtonian expansion ceases to be accurate. 

• We restrict ourselves to circular orbits by omitting 
the orbital eccentricity; this hypothesis is consid- 
ered rather realistic for the supermassive binary 
systems visible in the LISA band we consider. 

• The contributions of spins are negligible or they 
are oriented in such a way that no spin-induced 
precession of the orbital plane takes place; we take 
care of spin contributions only into the waveform 
phase. We choose this hypothesis in order to con- 
trol the complexity of the problem and focus on 
the comparison of the FWF versus RWF. 

• We approximate the waveform at the 2PN order, 
both in amplitude and in phase, considering up to 
six harmonics in the case of FWF. This simplifica- 
tion is motivated by computational reasons due to 
our limited computational resources and does not 
affect in any significant way the final results. 

• Out of the seventeen parameters on which the most 
general waveform depends on, the 2PN approxi- 
mation wc arc considering here depends only on 
eleven parameters: the luminosity distance, four 
angles defining the constant source position and 
orientation of the binary in the orbital plane, two 
mass parameters, two parameters related to the 
spin-orbit and spin-spin coupling and one arbitrary 
phase and time. 

We assume that the observations are carried out ac- 
cording to the following: 

• We consider sources at cosmological distances, as 
the event rate of massive black hole binary sys- 
tems in our local Universe is likely to be negligible 
and only taking into account the whole Universe it 
becomes of significant importance; moreover there 
is a great interest in carrying high redshift sur- 
veys; indeed, unless differently stated, we consider 
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fiducial sources to be at redshift z = 1, in a flat 
Universe described by the following cosmological 
parameters: Hq = 71 km s'-^ Mpc~^, flm = 0.27 
and f^A = 0.73; the corresponding luminosity dis- 
tance is therefore 



l + z 



dz' 



Ho Jo ^n„,{i + z'f + 



(5.1) 



All the parameters considered here are the observed 
ones; they differ from the values of the parameters 
as measured in the source rest frame according to 



/ 
t 

M 
M 
A* 



/ 



l + z 
{l + z)t 

(1 + z)M 

{l + z)M 

(1 + 



(5.2) 



• Unless differently stated, we consider that LISA 
observes the inspiral for a whole year before it 
reaches the ISCO. This corresponds to different 
frequency ranges depending on the harmonic. For 
every choice of source parameters these frequen- 
cies are computed in advanced and we also im- 
pose a low-frequency cut-off to the instrument at 
5 X 10-5 Hz. 

• The total noise that affects the observation is given 
by Eq. (|2.4p . that is, we take both the instrumental 
and confusion noise contributions. 

• We compute the expected mean square errors 
((AA*^)^)^/^ and the angular resolution of the in- 
strument, which we define as 



AQat = 27r^(Acos6i^)(A0^) - (A cos6'jvA0Ar)2 , 

(5.3) 

where {9^,4)^) are the polar angles in the solar 
system barycentre frame of the source location in 
the sky, with one and with both detectors, but we 
only report the results for the combined case. 

• The analysis is done in the frequency domain us- 
ing the stationary phase approximation: we first 
compute analytically the derivatives djh^'^\ where 
j = 1, . . . , 11, then compute numerically the Fisher 
matrix F^ and the variance-covariance matrix E-'*^; 
the integration and matrix inversion are performed 
using numerical routines of the gsl library. 

• To provide an overall picture of the instrument per- 
formances, we study not just a few cases in detailed 
but we also do extensive Monte-Carlo simulations, 
by varying the relevant source parameters, in par- 
ticular the position and orientation of the source, 
as it turns out to affect very significantly the pa- 
rameter measurements. 



B. The post-Newtonian inspiral waveforms 

In this section we derive explicit ready-to-use analyt- 
ical expressions for the signal measured at the LISA de- 
tector output for inspiral binary systems in circular orbit 
within the 2PN approximation, that could easily be ex- 
panded to include higher order terms. 

For our analysis, it is convenient to expand ft,^*' (t) as 
a summation of different multipole terms, which can be 
written schematically as 

i=i j=i ^ 



(u+jF+ 4- Uxj-Fx)cos(-$ + <^d) 



(5.4) 

where ipo is the LISA's Doppler phase, {it+_x , 'fi'-i-.x } 
contain an internal summation over all FN orders (see 
Eq. (|5.6p below), x is the FN expansion parameter 



{Mw) 



1/3 



(5.5) 



being w the orbital frequency and $ is given in Eq. (jB5p . 

In general, x)j an d can be written (at 

least, up to 2PN) as Eq. (|5.6p . Analyzing those terms, 
given explicitely in [H^l, one realizes that all of them 
have a common factor T(j) that only depends on the 
multipole j we are working with, regardless of the FN 
order considered 



E 



X u 



(n) 



T(j)E^" 



~(n) 



n=0 



(+,x),j 



n=0 

= T(j)^i(+,x),i 

E-"-(:!x). = T^wE-"-(+.x). 

n=0 

= T(j)u^(+,x)j • 



in) 



n=0 



(5.6) 



This notation is very convenient for computing analyt- 
ically the derivatives of the waveform. All (u, ?i))|"'' . 
can be found in Appendix IB 21 The factors T(j), for 
each multipole, are 



T(.? 


= 1) ^ 


Sm 


T(j 


= 2) ^ 


1 


T(j 


= 3) ^ 


Sm 

'JT 


T(j 


= 4) ^ 


s' 


T(j 


= 5) ^ 


n Sm 


T(j 


= 6) ^ 





(5.7) 
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where s = sint =| L x N | and c = cosi = — L • N. 
The source location in the sky N and the orbital angu- 
lar momentum L can be described by the polar angles 
(^jv,</'Ar) and {Ol,4>l) with respect to the solar system 
barycentrc frame. 

With all these considerations, the gravitational wave- 
form, given by Eq. (|5.4p . can be rewritten as follows 



"Pp.. 



where 



(5.8) 



A. 



I TO') I [{u+.,F++u^.,F^f 
+ {wxjF^ + w+jF+) 

' (Wx j-Fx + W+^jF+) 



tan 



(5.9) 



with respect to the different parameters, are 



where 



\ dSm ^ 
dM dM , 
dM d5m' 

\d5m J 

dM dSm 



(dM\ MM 
VdM) ^ 



(6.2) 



[m) ^ 



(6.3) 



dSm dSm 



dM dM 



f dMdn dM d^L \ ^Sm 

\ dM dSm dSm dM J ^ ' ' 



dM M 



dM 20fi 



(6.5) 



Note that T(j) cancels in Lpp,j expression, which prevents 
a lot of divergence problems in the numerical computa- 
tion of the Fisher matrix. Sec Appendix IB II for further 
details regarding the waveform in the frequency domain. 



dM 

dSm 



dfi 
dM 



3M Sm 
10m m 



4M2 



(6.6) 



(6.7) 



VI. RESULTS 

The signal model considered here depends on eleven 
independent parameters, and it is rather natural to con- 
sider the choice 

A = {cos8N,4>N,cos9L,4'LA'o^DL,tc,4>c, f3,(J, lnM,\nfi} 

(6.1) 

that we adopt in this study, in order to easy comparison 
with previous existing results. <c and (f>t are the time 
and phase at coalescence, and /3 and a arc the so called 
spin-orbit and spin-spin parameters. We usually set tc = 
(j)c = P = (J = in our analysis. 

The code we use to compute the Fisher matrix is an 
extension of a previous one developed by Alberto Vcc- 
chio. Our code allows for multiple choices: selection of 
PN order in amplitude and phase, number of harmon- 
ics, and also different sets of independent parameters. 
In particular, for the FWF it is more convenient to con- 
sider a different choice of the mass parameters to avoid 
divergences. Therefore, sometimes we consider the mass 
parameter combination {dm, Af }, where Sm ~ mi — TO2 
instead of {M,ij,}, although we always express the re- 
sults for comparison with respect to the later choice. The 
relation of the variance-covariance matrix components. 



Unfortunately, for the equal mass case, the Jacobian of 
the transformation between {dm, M} and {M, fi} is sin- 
gular, and also the Fisher matrix presents a coordinate 
singularity depending on the choice of mass parameters 
and the waveform model used. Since we are still inter- 
ested in the errors in In M and In ^ we convert them us- 
ing the same previous formulas. For the unequal masses, 
we find that computing the errors in {Sm, M} and then 
converting gives the same result as simply computing the 
errors in {M , yu} directly. Further details are discussed 
in Appendix [XI 

We have checked that our results agree with Vccchio's 
code for the RWF, and with Sintes and Vecchio at 
the 0.5PN-2PN order in amplitude and phase, respec- 
tively. We use numerical integration and matrix inver- 
sion routines from the gsl library that, for some partic- 
ular cases, we have checked against Mathematica. 

We work with the full 11 x 11 Fisher matrix for both 
RWF and FWF. There are parameter configurations for 
which cither the numerical integrations or the inversion 
of such a large matrix tend to fail (or the results do not 
have the desired accuracy). The reason is that the Fisher 
matrix is often ill conditioned. Because of that, we have 
limited our study up to a mass ratio m2/mi = 0.01. 
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10^ 10* 10** 10*' 10' 10* 

Total Mass (IV^yf^) 



FIG. 1: SNR versus total mass for the mass ratios of 1, 0.1 
and 0.01 for an integration time of one year. The sohd hnes 
correspond the FWF and dashed hnes to RWF. The sources 
are at redshift z = 1, corresponding to a luminosity distance 
of Dl = 6.64 Gpc, with fixed angles given by cos 6n = —0.6, 
<j)N = 1, cosOl ~ 0.2 and <j)L = 3. 



A. The impact of the FWF: General trends 

Given the extensive parameter space to be explored, 
we start by considering sources with a fixed location and 
orientation given by cos 6* AT = —0.6, 4>n = 1, cos9l = 0.2 
and (j)L — 3. This is the same choice as in [2^, and, as 
we will point out later, it corresponds to a case for which 
the SNR and the mass errors are similar to the typical 
average ones. Note that best and worst cases can span 
several orders of magnitude. For this set of angles, we 
consider sources at redshift z ~ I, corresponding to a 
luminosity distance of Dl = 6.64 Gpc and we study the 
effects on the SNR and parameter estimation using the 
RWF and the FWF as function of the binary total mass. 

In Fig. [T] we plot the SNR computed using the RWF 
and the FWF as a function of the total mass of the binary 
system. For systems whose total mass Af < 4 x 1O''M0 
the RWF over-estimates in general the SNR by a few 
percent. This was already pointed out in [2^. It is also 
a known fact that the Newtonian amplitude is about 7% 
higher than the 2.5PN order amplitude and also than 
the amplitude obtained by numerical simulations over 
the last few orbits before merger. For binaries with AI > 
4x K^Mq the second harmonic is no longer visible in the 
LISA band and higher harmonics, and therefore the use 
of FWF, play an important role extending the mass reach 
for supermassive black holes. The 'jumps' at high masses 
are due to the low-frequency cut-off to the instrument at 
5 X 10~^ Hz. In the case of the FWF we have limited 
our study to systems with M > 2 x IO^Mq due to our 
limited computational resources. 

For the same configuration, in Fig. [5] we represent the 



errors of the most relevant parameters as function of the 
total mass for different mass ratios. The Fisher matrix 
has been evaluated assuming the black hole spins to be 
zero, so that the spin-orbit and spin-spin parameters, 
P and cr, respectively, are equal to zero. In all cases, 
the errors are smaller for the FWF. For a total mass 
M < lO^Af0 the improvements are modest, except for 
the mass estimation for nearly equal masses, for which 
the errors in M and fx are of the same order; while these 
improvements are considerable for M > 5x 10^ Mq. This 
betterment is not due to an increase of SNR, but to 
the higher harmonics that contribute to disentangle the 
source parameters. For the equal mass case, the errors in 
/X improve up to 3 orders of magnitude at M = IO^AIq. 
In general the measurements of the masses improves by 
more than an order of magnitude for M > IQ^AIq. We 
find interesting to add the case TO2/mi = 0.9 obtained 
by the FWF (dotted dark line, in the lower right panel) 
for comparison. Using the FWF the luminosity distance 
can be measured to better than 10% up to 3 — 4 x IO^Mq 
for this particular source location, while using RWF the 
errors are several orders of magnitude larger for these 
large masses. Similarly, the error box in the sky improves 
significantly at 1O''M0, although this sky position, as we 
will see below, was not one of the most favorable ones. 

Using the FWF and a mass ratio of 1712/1711 = 0.01, 
there were several cases in which the errors we obtained 
did not have the desired accuracy. Those cases have been 
marked in Fig. [2] by replacing the solid line by a dotted 
one. 

LISA will also be able to measure gravitational waves 
from massive black hole coalescences to large redshift, 
making possible to study the merger history of black 
holes. Therefore, it is interesting to extend the anal- 
ysis to higher redshifts. Figure [3] shows the luminos- 
ity distance as a function of redshift for a flat Uni- 
verse described by the cosmological parameters: Hq = 
71 km s-i Mpc"\ = 0.27 and = 0.73. In Fig- 
ures m and [5] we plot the SNR versus total mass for 
redshifts z ~ 1, 10 and 20. The results we obtain are 
the expected ones since, modulo over-all amplitude, the 
gravitational waves that we measure from a binary with 
masses {A4 , 11} at redshift z are those of a local system 
with masses {(1 -|- z)A4, (1 -I- z)fi}. The SNR decreases 
with z, not only due to the distance, but also because 
the frequency of the signal is redshifted and the total 
effective noise for LISA is higher at lower frequencies. 

In Figure [6] we plot the distance measurement er- 
ror, angular resolution, and mass measurement errors 
for LISA observations of the final year of supermassive 
black hole inspirals. The fiducial sources are at z =1, 
10 and 20. The waveform considered are the FWF at 
2PN order and the RWF with cos6'Ar = -0.6, 0^ = 1, 
cos6'l = 0.2, = 3 and f3 — a ^ 0. For simplicity 
we display the curves corresponding to the equal mass 
case, but, as expected, other mass ratios follow the same 
trend. All curves drift with the redshift, but qualita- 
tively the impact of FWF versus RWF is the same as for 
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FIG. 2: Errors versus total mass for the same case as Fig. [T] The sohd hnes correspond to FWF and dashed hnes to RWF. 
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FIG. 3: The luminosity distance as a function of redshift 
for a flat Universe described by the cosmological parameters: 
Ho = 71 km s"^ Mpc~\ Sl„ = 0.27 and f^A = 0.73. 
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FIG. 4: SNR versus total mass for various values of the source 
redshift. We consider equal mass binary systems with fixed 
angles given by cos 9n ~ —0.6, 0jv ~ 1, cos 9l = 0.2 and 
(j)L = 3. The solid lines correspond the FWF and dashed 
lines to RWF. 



z = 1, but for different mass ranges. 



B. Exploring the parameter space 

As we mentioned in the previous section, the error 
measurements are very sensitive - they vary by orders-of- 
magnitude - to the true value of the source parameters; 
in order to give meaningful results, one therefore is forced 
to explore a large parameter space. We do this (i) by 
considering sources on an isotropic grid in the sky with 
a fixed orientation, and (ii) by extensive Monte-Carlo 
simulations for all possible location and orientation of 



z = 1 
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FIG. 5: The same as Fig. |4]for a mass ratio m,ilm\ — 0.01. 



the source with respect to LISA. 

Figure [7] is a sample of a sky map of SNR for LISA ob- 
servations of the final year of inspirals using FWF. The 
sources considered correspond to toi = m2 = IO^Mq 
at redshift z = 1, with orientation angles cos 9l = 0.2 
and (j)L — 3. For all sources we assume the location 
of LISA at the time of coalescence is 4'lisa = 0. The 
SNR over the entire sky covers a range ~ 100 — 900. The 
SNR is higher for sources located orthogonal to the plane 
of LISA at the time of coalescence, which corresponds 
to ((/)jv = 0°;eN = 120°) and {(j^N = 180°; 61^ = 60°), 
or in general (0^ ^ 4>lisa;0±_ = 9 lis A + 30°) and 
(<^_L = cl)LisA + 180°; = Olisa - 30°), since the LISA 
constellation is inclined at an angle of 60° with respect to 
the ecliptic. The reason is that most of the SNR is accu- 
mulated in the last days before merger. The reader can 
notice also the nearly symmetric [{9, (p) (tt — 6*, + tt)] 
profiles of figures [TlfTUl The small asymmetry shows the 
relative importance of the Doppler phase modulation. 

Figure [8] shows the angular resolution, and the error in 
the measurements of the luminosity distance and masses 
over the entire sky for the same case as Fig. [T] To sum- 
marize, the parameters of two 10'' black holes spiral- 
ing toward the final merger at z = 1 can be measured 
very accurately depending on the sky location: Arjjv up 
to ~ 10~^ srad, the luminosity distance to better than 
1%, and the masses between 1% and 10%. Although 
it is true that for a given source, the higher SNR the 
better the parameter estimation, e.g., by changing the 
distance, this cannot be generalized comparing the SNR 
at different locations in the sky. It is not just a matter 
of SNR but long observation times that contribute to 
disentangle and improve the parameter estimation. 

In Fig. [9] we show the gain in angular resolution com- 
paring the results obtained using FWF versus RWF for 
the case mi = m2 = IO^Mq as in Fig. [S] The benefit of 
using FWF is clear. For all those sky locations in which 
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Total Mass (Mg^i^) Total Mass (l^y^^) 

FIG. 6: Distance measurement error, angular resolution, and mass measurement errors for LISA observations of the final 
year of equal mass supermassive black hole inspirals. The fiducial sources are at z =1, 10 and 20. The waveform considered 
are the FWF at 2PN order (solid lines) and the RWF (dashed lines) with cos 6n = —0.6, </)jv = 1, cos 9l = 0.2, (f>L = S and 
/3 = (7 = 0. 
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FIG. 7: Sky map of SNR for LISA observations of the final year of inspirals using FWF. The sources considered correspond 
to mi = 7712 = W^Mq, at redshift z — 1, with orientation angles cosOl = 0.2 and cj>L ~ 3. For all sources we assume the 
location of LISA at the time of coalescence is 4>lisa = 0. 
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FIG. 8: Sky maps for the angular resolution, distance measurement error and mass measurement errors for LISA observations 
of the final year of inspirals using FWF. As in Fig.[7]these correspond to mi = m2 — 10^ Mq at redshift z — 1, with orientation 
angles cos Ol = 0.2, 0_l = 3 and f3 — a = 0. For all sources we assume the location of LISA at the time of coalescence is 
(pLiSA = 0. The dark blue dot corresponds to cos 9n ~ —0.6, <j!>jv = 1, mentioned in the text. 




FIG. 9: Gain in angular resolution comparing the results obtained using FWF versus RWF for the equal mass case mi = 
m2 = IO^Mq with the same assumptions as in Fig. (8] The dark blue dot corresponds to cos 9n = —0.6, 4>n = 1- 
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FIG. 10: Sky map of the gain in angular resolution for LISA observations of the final year of inspirals using FWF versus 
RWF corresponding to mi — IO^Mq, m2 = 10^ Mq and z = 1. We assume all source have the same orientation cos 6l ~ 0.2, 
(l)L ~ 3, zero spins fi — a — and that LISA is at cj>LiSA = at the time of coalescence. The dark blue dot corresponds to 
cosSjv = —0.6, (/)]v = 1. 



we obtain the best angular resolution, AfJjv ^ 10~^ srad 
with the FWF, the corresponding gain is up to 3 orders 
of magnitude. 

Figure [TO] represents also the improvement in angu- 
lar resolution for the unequal mass case mi = IO^A/q, 
m2 = 10® M0. In this case the optimal gain is of ~ 2 
orders of magnitude, and those correspond to the same 
sky locations as in the equal mass case analyzed before. 

In order to cover completely the full parameter space, 
we proceed to perform Monte-Carlo simulations accord- 
ing to the following: we consider an ensemble of fiducial 
sources all at redshift z = 1 (which sets the luminosity 
distance Dl = 6.64 Gpc), with zero spins /3 = = 
and we select the value of the masses mi and m2. For 
each set of mass parameters we select randomly the four 
geometrical angles (^at, (pNi Ol and (Pl) from an uni- 
form distribution in cosOn, (t>N, cosOl and 0l, and as 
far as for the other two parameters, we chose them as 
tc = 4>c = 0. The Monte-Carlo is done on a 1000 dif- 
ferent sets of angles. We present the results in terms of 
probability distributions. 

We have studied six different pairs of masses: (a) mi = 
IO^Mq; ma = IO^A/q, (b) mi = IO^Mq; ma = 10® AT©, 
(c) mi = IO'^ATq; ma = IQ^Mq, (d) mi = lO^Afo; 
ma = 10*^ Mq, (e) mi = IO^Mq; ma = IO'^Mq, and 
(f) mi = IO^Mq; ma = 10'"^ Mq. Figure [ll] and table |T] 
summarize the results. 

The key result, is that using the FWF the errors are 
smaller than with the RWF. There are big improvements 
for the IO^ATq - lO^Mg, and 10^ Af© - lO^Af© in angu- 
lar resolution and distance measurement: the angular 
resolution improves in average 25 and 7.3, respectively; 
and the luminosity distance by factors of 62 and 6.7, 
respectively. One should also notice that, in those two 
cases, those parameters were poorly determined using 



only the RWF. For the other sets of masses the aver- 
aged improvement in angular resolution and luminosity 
distance are more moderate, between 1.3 — 3 for ASIat, 
and 1.7 — 2.7 for ADi/D^. In all cases the masses are 
determined much more accurate, even by several orders 
of magnitude in the case of /i, using the FWF. For the 
equal mass cases, the errors in Ai and /i are of the same 
order using the FWF. 

Because of the different harmonics, the FWF has a 
much greater richness than RWF that clearly improves 
the parameter estimation. It is worth mentioning that 
similar level of improvements were obtained in [2^ where 
the waveform considered was only at the 0.5PN-2PN or- 
der in amplitude and phase, respectively, i.e., adding the 
first and third harmonics. This suggests that the im- 
provement in parameter estimation is mainly due to the 
inclusion of the third harmonic, which also increases the 
mass reach of LISA. The importance of the different PN 
orders has been discussed in detailed recently by Arun 
et al. in [H, [l^] for so me p articular cases, and also for 
ground-based detectors [3g |. 

In Figure [T^] we compare the distribution of the mea- 
surement errors for a waveform containing only the sec- 
ond and third harmonics, keeping both amplitude and 
phase at the 2PN order, with the FWF and the RWF, 
for the pair of masses mi = l{f Mq and ma ~ 10® Af©. 
This figure shows how a substantial improvement in pa- 
rameter estimation is obtained by adding the third har- 
monic. In this case, the mean and standard deviation 
of the logarithm of the measurement errors are the fol- 



lowing: logiQ Af2 



'TV 



-2.33 ± 1.06 srad, lognjAr^L 



-2.05 ± 0.46 srad, log^^ I^Dl/Dl 



-1.22 ± 0.27, 
-2.44 ±0.17 

and logiQ A/i//i = — 1.57± 0.36. These results are very 
close to those obtained for the FWF as can be seen from 



logio Ate = 2.96 ± 0.19 s, logn, ATW /M 
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FIG. 11: The probability distributions of SNR and measurement errors for observations of the final year of supermassive black 
hole binaries at z = 1 by LISA. For each pair of masses, the histograms show the result of two Monte Carlo simulations, where 
1000 sources have been randomly located and oriented in the sky. The plots compare SNR and the errors for the FWF (in 
clear green) and the RWF (in dark red). 
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TABLE I: Characterization of the probability distributions of SNR and measurement errors of Fig. 1111 For each pair of masses 
and waveform model used, the mean and standard deviation of the SNR is given, as well as the mean and standard deviation 
of the logarithm of the measurement errors. The gain factors are computed as {SNR)fwF/{SNR)nv/F and lo'<^>i^«F~^^>Fwf)^ 
for the SNR and measurement errors, respectively. 
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FIG. 13: The progressive accumulation ol SNR as a lunction 
of a look-back time; The observations refer to the final year of 
the inspiral of supermassive black hole as recorded by LISA 
for fiducial sources at redshift z = 1, with cosOn = —0.6, 
(j>N = 1, cosOl ~ 0.2 and (pL ~ 3. The solid lines correspond 
the FWF and dashed hues to RWF. 



FIG. 12: Comparison of measurement errors for the FWF 
(in clear green), the RWF (in dark red) and a waveform con- 
taining only the second and third harmonics (in very dark 
blue), corresponding to LISA observations of the final year 
of inspirals for mi — IO^Mq, m2 — 10® and z = 1. 



table |T] and are also in agreement with the level of im- 
provement found in [2^ when considering only the 0.5 
PN order in amplitude. 



C. Pre-merger localization 

From an astronomical point of view, one of the most 
attractive features is the possibility that LISA might 
have enough angular resolution to locate the galaxy or 
galaxy cluster where the coalescence of a massive black 
hole takes place and therefore identify potential electro- 
magnetic counterparts. The angular resolution is de- 
duced primarily from the detector's motion around the 
Sun, so one expects that the uncertainty in the angu- 
lar resolution will not change so much during the last 
days before merger. Therefore, as discussed in [l3|, we 
are interested in analyzing the time-dependence of the 
angular resolution and SNR, as a function of some look- 
back time tisco ^ tf prior to coalescence, and measure 
the importance of the FWF versus RWF. 

For fiducial sources at z — I and a given sky lo- 
cation and orientation. Figures [13] and [M] confirm the 
importance of the FWF, in particular for those sys- 
tems with a higher total mass: IO'^Mq — 10'' and 
lO^M^o — lO^Mrri, while there is not much difference for 




FIG. 14: Time-dependence angular resolution as a function 
of a look-back time for the same sources as in Fig. [TS] The 
solid lines correspond the FWF and dashed lines to RWF. 



the lO^M© — IO^A/q one. For equal masses, only the 
even multipoles (2, 4, 6) contribute to the FWF, while 
for the unequal masses there are contribution from all 
the six harmonics. The 'jumps' in the progressive accu- 
mulation of SNR correspond to those times in which a 
new higher harmonic enters into the LISA band and it 
is related to the lower frequency cut-off at 5 x 10~^ Hz 
we have imposed. For example, in the 10'' Mq — lO^M© 
case the contribution of the 4th harmonic becomes rel- 
evant around 10 days before coalescence while the 2nd 
harmonic rapidly increases the SNR 2 days before coa- 



16 



lescence. For the unequal mass case 10'' Af© — IO^Mq 
we clearly see the contributions of the 2nd, 3rd and 4th 
harmonics. 

Using the FWF, not only the uncertainty in sky loca- 
tion decreases but also allows earlier warnings. 



VII. SUMMARY AND OUTLOOK 

We have considered LISA observations of supcrmas- 
sive black hole systems in the final stage of inspiral. We 
have restricted our analysis to systems in circular orbit 
with negligible spins, modeling the radiation at the full 
2PN order, and we have compared with the restricted- 
2PN. With both waveform models wc have determined 
the mean-square errors associated with the parameter 
measurements of black hole binaries in the mass range 
IO^Mq — IO^Mq, for equal and unequal mass cases, for 
a wide range of source locations and orientations. 

The conclusions of this work are particularly impor- 
tant with regard to the astrophysical reach of future 
LISA measurements. Our analysis clearly shows that 
modeling the inspiral with the full post-Newtonian wave- 
forms, as compared to the restricted-PN ones, not only 
extends the reach to higher mass systems up to 10^ Mq, 
as previously shown in [26[ , but also improves in general 
the parameter estimation, and allows for early warnings 
for systems with a high total mass. There are remark- 
able improvements in angular resolution and distance 
measurement for systems with a total mass higher than 
5 X IO^Mq, as well as a large improvement in the mass 
determination. For A^/^, the improvement is more than 
three orders of magnitude in the case of equal masses. 

For binary systems of lO'^M© - lO^Mg, and IO^Mq - 
IO^Mq at redshift z = 1, the angular resolution improves 
in average 25 and 7.3, respectively; and the luminosity 
distance by factors of 62 and 6.7, respectively. Moreover, 
for the equal mass case IO^Mq — IO'^Mq, for all those 
sky locations in which we obtained the best angular res- 
olution, Ar^AT ~ 10~^ srad with the FWF, the gain in 
AfZ^v is up to 3 orders of magnitude. These results are 
in agreement with those recently found in p7| . 

These improvements are related to the fact that the 
FWF has a much greater richness than the RWF, due to 
the presence of the higher harmonics and, in particular, 
the main contribution to this improvement is associated 
to the third harmonic. 

There are a number of issues that influence these ob- 
servations: the instrumental lower-frequency cut-off we 
have imposed and the confusion noise model we have 
used are two of them. It would also be very interesting 
to revise these conclusions, by considering black holes 
with large spins and precession, since it is known that 
the presence of spins reduces the errors with which the 
source parameters are measured [T9| . 

Another issue is the fact that the largest improvement 
happens for systems with a total mass of 5 x lO^Af© and 
higher, and although, we are still in a regime of large 
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FIG. 15: Mass measurements using the FWF for binary sys- 
tems with a total mass M — 10^ Mq. The solid lines cor- 
respond to compute the errors in {5m, M} and converting 
them, by mean of equations (|6.2|l - (|6.8|) . into {M,^} , and 
crosses correspond to compute the errors in {A4, n} directly. 



SNR, one could question the validity of the Fisher matrix 
approach. Other investigations are currently underway 
using alternative methods [6l| . 



APPENDIX A: MASS PARAMETERS 
TRANSFORMATION 



As discussed earlier in Sec. I VII it is traditional to pa- 
rameterize the masses using In Ai and In /i because of 
their appearance in the waveform phase. However the 
higher order PN amplitude terms depend on the mass 
difference Sm and total mass M. It is a matter of choice 
to work with one or another set of parameters. Errors in 
{M,fj.} can be computed easily, in principle, given the 
errors in {Sm, AI } by mean of equations (|6.2p - (|6.8[) . Un- 
fortunately the Jacobian of the transformation between 
{A^,/i} and {5m, M} is singular when dm ~ leading 
to problems in evaluating the Fisher matrix. 

For unequal masses, we find that computing the errors 
in {6m, M} and then converting gives the same result 
as simply computing the errors in {M , fi} directly, and 
this result is independent of the waveforms used FWF or 
RWF, as it is shown in figures [TSl and [TBI However, de- 
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pending of our choice of mass parameters and waveform 
model we use, there appear divergences when evaluat- 
ing the Fisher matrix for equal masses. In particular, 
for the RWF we do not trust the {Sm, M} parameter- 
ization, and for the FWF we do not trust the {Ai,^}. 
That is, wc trust the solid lines in Fig. [15] and [161 but 
not the crosses for equal masses, that either do not coin- 
cide with the solid lines for the FWF or could not even 
be computed using the RWF. Because the solid lines did 
not present any misbehavior nor discontinuities when ap- 
proaching the equal mass case, we decided to use (inde- 
pendently of the masses) the {Ai , fj,} parameterization 
for the RWF and the {Sm, M} for the FWF. 

Another aspect we want to study in more detail, is 
the fact that using the FWF, for equal masses, A M/Ai 
and A fi/^ become identically the same, but this is not 
true for the RWF, nor for unequal masses. This is a 
consequence of setting Sm = in Eqs. (|6.2p - (|6.8p . In 
Figurc[T7]we show that this is not an artifact of our mass 
transformation due to the Jacobian being singularity at 
that point but a fact. 



APPENDIX B: SOME EXPRESSIONS RELATED 
WITH THE 2PN INSPIRAL WAVEFORM 



FIG. 16: Mass measurements using the RWF for binary 
systems with a total mass M — lO^M©. The solid lines 
correspond to compute the errors in {Ai, fJ.} directly, and 
crosses correspond to compute the errors in {5m, M} first 
and converting them into {Ai,j-i}. 
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FIG. 17: Study of the mass parameter measurements in the 
limit 6m using the FWF with the {Sm, M} mass param- 
eterization for systems with a total mass (top) M — 10^ Mq 
and (bottom) M = IO^Mq. In both panels the solid lines 
correspond to A the dashed lines to A fi/jj, and the 

circle indicates the equal mass case. 



1. Time, frequency and phase evolution 

A coalescing binary system evolves by loosing energy 
and angular momentum L through emission of gravi- 
tational waves of increasing frequency and amplitude. 
Working with the post-Newtonian approximation, tak- 
ing into account possible spin motions of each object. Si 
and S2, the signal frequency of the second harmonic, 
F = 2forb, evolves, up to 2PN order, according to 
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where x is the PN expansion parameter defined in 
Eq. (|5.5[) . and f3 and a are the so called spin- orbit and 
spin-spin parameters, respectively 
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Integrating Eq. (jBl[) . one can derive the time evolution 
of the gravitational radiation 
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and the phase evolution of the gravitational waveform. 
For the second harmonic this is 
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The gravitational waveform can be computed in the 
frequency domain using the stationary phase approxi- 
mation 
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{dF/dt) ' can be written to 2PN order as 
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where hj = -^2Mr]-—x^Aj. 



Using Eq. dm]), 
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where coefficients fc„ are defined as follows 
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Contributions to w^^. 
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and 77)(+ x).j up to 2PN 
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In Sec. IV Bl we have seen a general form to expand 
the GW amplitude, of a particular multipole, as a sum- 
mation of its different PN contributions. The analytical 
expression of all the terms appearing in Eq. (|5.6p can 
be obtained from [HO]- Here we explicitely give all the 
non- vanishing terms uS"'„s and wf"''„s up to ti = 4. 

Contributions to u+j 
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